o 
o 




NT-UW 05-13 
CALT 68-2579 



Ginsparg-Wilson Pions Scattering in a Sea of Staggered Quarks 

Jiunn-Wei Chen/'Q Donal 0'Connell,^'lll Ruth Van de Water, ^'P and Andre Walker-Loud^ 'P 

^Department of Physics, National Taiwan University, Taipei 10617, Taiwan 

\q'' ^California Institute of Technology, Pasadena, CA 91125, USA 

O ■ ^Department of Physics, University of Washington, 
(N ■ 

^ ■ Box 351560, Seattle, WA 98195-1560, USA 

<^ ! Abstract 

' We calculate isospin 2 pion-pion scattering in chiral perturbation theory for a partially quenched, 

■ mixed action theory with Ginsparg-Wilson valence quarks and staggered sea quarks. We point 

> : 

Tlj" . out that for some scattering channels, the power-law volume dependence of two pion states in 

\ nonunitary theories such as partially quenched or mixed action QCD is identical to that of QCD. 



Thus one can extract infinite volume scattering parameters from mixed action simulations. We 
then determine the scattering length for both 2 and 2+1 sea quarks in the isospin limit. The 
scattering length, when expressed in terms of the pion mass and the decay constant measured on 
the lattice, has no contributions from mixed valence-sea mesons, thus it does not depend upon 
the parameter, Cmix; that appears in the chiral Lagrangian of the mixed theory. In addition, 
^ , the contributions which nominally arise from operators appearing in the mixed action 0{a^mq) 

\ Lagrangian exactly cancel when the scattering length is written in this form. This is in contrast 

to the scattering length expressed in terms of the bare parameters of the chiral Lagrangian, which 
explicitly exhibits all the sicknesses and lattice spacing dependence allowed by a partially quenched 
mixed action theory. These results hold for both 2 and 2+1 fiavors of sea quarks. 
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I. INTRODUCTION 



Lattice QCD can, in principle, be used to calculate precisely low energy quantities in- 
cluding hadron masses, decay constants, and form factors. In practice, however, limited 
computing resources make it currently impossible to calculate processes with dynamical 
quark masses as light as those in the real world. Thus one performs simulations with quark 
masses that are as light as possible and then extrapolates the lattice calculations to the 
physical values using expressions calculated in chiral perturbation theory (xPT). This, of 
course, relies on the assumption that the quark masses are light enough that one is in the 
chiral regime and can trust xPT to be a good effective theory of QCD 



can at present reach significantly lighter 



Lattice simulations with staggered fermions 
quark masses than other fermion discretizations and have proven extremely successful in 
accurately reproducing experimentally measurable quantities Staggered fermions, 

however, have the disadvantage that each quark flavor comes in four tastes. Because these 
species are degenerate in the continuum, one can formally remove them by taking the fourth 
root of the quark determinant. In practice, however, the fourth root must be taken before 
the continuum limit; thus it is an open theoretical question whether or not this fourth- 
rooted theory becomes QCD in the continuum limit. ^ Even if one assumes the validity of 
the fourth-root trick, which we do in the rest of this paper, staggered fermions have other 
drawbacks. On the lattice, the four tastes of each quark flavor are no longer degenerate, and 
this taste symmetry breaking is numerically significant in current simulations p . Thus one 
must use staggered chiral perturbation theory (S^PT), which accounts for taste-breaking 
discretization effects, to extrapolate correctly staggered lattice calculations to the continuum 
0, y, ^, 0- Fits of SxPT expressions for meson masses and decay constants have been 
remarkably successful. Nevertheless, the large number of operators in the next-to-leading 
order (NLO) staggered chiral Lagrangian Q arrd the compheated form of the Uaorr B- 
parameter in S^PT [11] both show that S^PT expressions for many physical quantities will 
contain a daunting number of undetermined fit parameters. Another practical hindrance to 
the use of staggered fermions as valence quarks is the construction of lattice interpolating 
fields. Although the construction of a staggered interpolating field is straightforward for 



^ See Ref. |y| for a recent review of staggered fermions and the fourth-root trick. 
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mesons since they are spin objects [12, ll3|, this is not in general the case for vector 



mesons, baryons or multi-hadron states since the lattice rotation operators mix the spin, 



15 



angular momentum and taste of a given interpolating field jl^ . 

The use of Ginsparg- Wilson (GW) fermions j3| evades both the practical and theoretical 
issues associated with staggered fermions. Because GW fermions are tasteless, one can 
simply construct interpolating operators with the right quantum numbers for the desired 
meson or barvon. Moreover, massless GW fermions possess an exact chiral symmetry on 
the lattice [l8| which protects expressions in xPT from becoming unwieldy.^ Unfortunately, 
simulations with dynamical GW quarks are approximately 10 to 100 times slower than 
those with staggered quarks j^l and thus are not presently practical for realizing light 
quark masses. 

A practical compromise is therefore the use of GW valence quarks and staggered sea 
quarks. This so-called "mixed action" theory is particularly appealing because the MILC im- 
proved staggered field configurations are publicly available. Thus one only needs to calculate 
correlation functions on top of these background configurations, making the numerical cost 
comparable to that of quenched GW simulations. Several lattice calculations u sing domain- 



wall or overlap valence quarks with the MILC configurations are underway 



including a determination of the isospin 2 (/ = 2) vrvr scattering leng th [29|. Although this 
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27, m 



is not the first / = 2 tttt scattering lattice simulation 30], |3l|, |3^ |33|, , it is the only one 
with pions light enough to be in the chiral regime P, J Its precision is limited, however, 
without the appropriate mixed action xPT expression for use in continuum and chiral ex- 
trapolation of the lattice data. With this motivation we calculate the / = 2 vrvr scattering 
length in chiral perturbation theory for a mixed action theory with GW valence quarks and 
staggered sea quarks. 
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Mixed action chiral perturbation theory (MAvPT) was first introduced in Refs. 

nn 

1361 . 13 7| and was extended to include GW valence quarks on staggered sea quarks for both 



mesons and baryons in Refs. 



anc 

n 



continuum, infinite-volume xPT |40L 



41^. 



3£ |, respectively, tttt scattering is well understood in 
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45L I4(t| , and is the simplest two-hadron 



process that one can study numerically with LQCD. We extend the NLO xPT calculations of 



^ In practice, the degree of chiral symmetry is hmited by how well the domain- wall fermion |19l 1201 121| is 
realized or the overlap operator [2^|23,|2J| is approximated. 
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Refs 



4^ to MAy 



QCD (PQQCD) 
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T. A mixed action simulation necessarily involves partially quenched 
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in which the valence and sea quarks are treated 
differently. Consequently, we provide the PQxPT tttt scattering amplitude by taking an 
appropriate limit of our MA^PT expressions. In all of our computations, we work in the 
isospin limit both in the sea and valence sectors. 



This article is organized as follows. We first comment on the determination of infinite 
volume scattering parameters from lattice simulations in Section |nj focusing on the applica- 
bility of Liischer's method [ssl . 5^ to mixed action lattice simulations. We then review mixed 



action LQCD and MA^PT in Section UTTl In Section HVl we calculate the / = 2 vrvr scattering 
amplitude in MA^PT, first by reviewing vrvr scattering in continuum SU{2) xPT and then 
by extending to partially quenched mixed action theories with Nj = 2 and Nf = 2 + 1 
sea quarks. We discuss the role of the double poles in this process 0] and parameterize 
the partial quenching effects in a particularly useful way for taking various interesting and 
important limits. Next, in section |3 we present results for the pion scattering length in 
both 2 and 2 + 1 flavor MA^PT. These expressions show that it is advantageous to fit to 
partially quenched lattice data using the lattice pion mass and pion decay constant mea- 
sured on the lattice rather than the LO parameters in the chiral Lagrangian. We also give 
expressions for the corresponding continuum PQxPT scattering amplitudes, which do not 
already appear in the literature. Finally, in Section IVII we briefly discuss how to use our 
MAxPT formulae to determine the physical scattering length in QCD from mixed action 
lattice data and conclude. 



II. DETERMINATION OF SCATTERING PARAMETERS FROM MIXED AC- 
TION LATTICE SIMULATIONS 

Lattice QCD calculations are performed in Euclidean spacetime, thereby precluding the 
extraction of S-matrix elements from infinite volume Liischer, however, developed a 

method to extract the scattering phase shifts of two particle scattering states in quantum field 
theory by studying the volume dependence of two-point correlation functions in Euclidean 
spacetime [s^l. In particular, for two particles of equal mass m in an s-wave state with 
zero total 3-momentum in a finite volume, the difference between the energy of the two 
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particles and twice their rest mass is related to the s-wave scattering length:^ 



AEa : 



In the above expression, oq is the scattering length (not to be confused with the lattice 
spacing, a), L is the length of one side of the spatially symmetric lattice, and Ci and C2 are 
known geometric coefficients.^ Thus, even though one cannot directly calculate scattering 
amplitudes with lattice simulations, Eq. (^, which we will refer to as Liischer's formula, 
allows one to determine the infinite volume scattering length. One can then use the expres- 
sion for the scattering length computed in infinite volume xPT to extrapolate the lattice 
data to the physical quark masses. 

Because Liischer's method requires the extraction of energy levels, it relies upon the exis- 
tence of a Hamiltonian for the theory being studied. This has not been demonstrated (and 
is likely false) for partially quenched and mixed action QCD, both of which are nonunitary. 
Nevertheless, one can calculate the ratio of the two-pion correlator to the square of the 
single-pion correlator in lattice simulations of these theories and extract the coefficient of 
the term which is linear in time, which becomes the energy shift in the QCD (and contin- 
uum) limit. We claim that in certain scattering channels, despite the inherent sicknesses 
of partially quenched and mixed action QCD, this quantity is still related to the infinite 
volume scattering length via Eq. ((T)), i.e. the volume dependence is identical to Eq. ((T)) up 
to exponentially suppressed corrections.^ This is what we mean by "Liischer's method" for 
nonunitary theories. We will expand upon this point in the following paragraphs. 

It is well known that Liischer's formula does not hold for many scattering channels in 
quenched theories because unitarity-violating diagrams give rise to enhanced finite volume 
effects 13]. For certain scattering channels, however, quenched xPT calculations in finite 



volume s 
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formula 
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61| 



. Chiral perturbation theory calculations additionally show that the 



same sicknesses that generate enhanced finite volume effects in quenched QCD also do so in 



partially quenched and mixed action theories 



Q H H H 0, 



It then follows 



^ Here we use the "particle physics" definition of the scattering length which is opposite in sign to the 
"nuclear physics" definition. 

^ This ex pres sion generalizes to scattering parameters of higher partial waves and non-stationary parti- 
cles [silsiElEi. 

^ Here, and in the following discussion, we restrict ourselves to a perturbative analysis. 
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that if a given scattering channel has the same volume dependence as Eq. in quenched 
QCD, the corresponding partially quenched (and mixed action) two-particle process will also 
obey Eq. Correspondingly, scattering channels which have enhanced volume dependence 
in quenched QCD also have enhanced volume dependence in partially quenched and mixed 
action theories. We now proceed to discuss in some detail why Liischer's formula does or 
does not hold for various 2^2 scattering channels. 

Finite volume effects in lattice simulations come from the ability of particles to propa- 
gate over long distances and feel the finite extent of the box through boundary conditions. 
Generically, they are proportional either to inverse powers of L or to exp(-mL), but Liischer's 
formula neglects exponentially suppressed corrections. Calculations of scattering processes 



in effective field theories at finite vo 
s-channel diagrams jsO 



61 
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ume show that the power-law corrections only arise from 
56^ . This is because all of the intermediate particles can 
go on-shell simultaneously, and thus are most sensitive to boundary effects. Consequently, 
when there are no unitarity-violating effects in the s-channel diagrams for a particular scat- 
tering process, the volume dependence will be identical to Eq. ((H), up to exponential correc- 
tions. Unitarity-violating hairpin propagators in s-channel diagrams, however, give rise to 
enhanced volume corrections because they contain double poles which are more sensitive to 



59|.^ Thus all violations of Liischer's formula come from on-shell hairpins 



boundary effects 
in the s-channel. 

Let us now consider / = 2 tttt scattering in the mixed action theory. All intermediate 
states must have isospin 2 and s > 4m^. If one cuts an arbitrary graph connecting the 
incoming and outgoing pions, there is only enough energy for two of the internal pions to 
be on-shell, and, by conservation of isospin, they must be valence vr'^^'s.'' Thus no hairpin 
diagrams ever go on-shell in the s-channel, and the structure of the integrals which contribute 
to the power-law volume dependence in the partially quenched and mixed action theories 



^ We note that, while the enhanced volume corrections in quenched QCD invalidate the extraction of scat- 
tering parameters from certain scattering channels, e.g. 1 = Isolls if. this is not the case in principle for 
partially quenched QCD, since QCD is a subset of the theory. Because the enhanced volume contribu- 
tions must vanish in the QCD limit, they provide a "handle" on the enhanced volume terms. In practice, 
however, these enhanced volume terms may dominate the correlation function, making the extraction of 
the desired (non-enhanced) volume dependence impractical. 

^ We restrict the incoming pions to be below the inelastic threshold; this is necessary for the validity of 
Liischer's formula even in QCD. 
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is identical to that in continuum xPT. This insures that Liischer's formula is correctly 
reproduced to all orders in 1/L with the correct ratios between coefficients of the various 
terms. Moreover, this holds to all orders in xPT, PQxPT, MAxPT, and even quenched 
xPT. The sicknesses of the partially quenched and mixed action theories only alter the 
exponential volume dependence of the 1 = 2 scattering amplitude.* This is in contrast 
to the J = vrvr amplitude, which suffers from enhanced volume corrections away from 
the QCD limit. In general, the argument which protects Liischer's formula from enhanced 
power-like volume corrections holds for all "maximally-stretched" states at threshold in the 
meson sector, i.e. those with the maximal values of all conserved quantum numbers; other 
examples include K'^K'^ and K^tt'^ scattering. We expect that a similar argument will hold 
for certain scattering channels in the baryon sector. 

Therefore the s-wave / = 2 vrvr scattering length can be extracted from mixed action 
lattice simulations using Liischer's formula and then extrapolated to the physical quark 
masses and to the continuum using the infinite volume MA^PT expression for the scattering 
length.^ 



III. MIXED ACTION LAGRANGIAN AND PARTIAL QUENCHING 

Mixed action theories use different discretization techniques in the valence and sea sectors 
and are therefore a natural extension of partially quenched theories. We consider a theory 
with Nf staggered sea quarks and Ny valence quarks (with corresponding ghost quarks) 



nee qua ] 

HQ. 



which satisfy the Ginsparg- Wilson relation [17|,|1^. In particular we are interested in theories 
with two light dynamical quarks {Nf = 2) and with three dynamical quarks where the 
two light quarks are degenerate (commonly referred to as Nf = 2 + 1). To construct the 
continuum effective Lagrangian which includes lattice artifacts one follows the two step 
procedure outlined in Ref. |68|. First one constructs the Symanzik continuum effective 
Lagrangian at the quark level IgqI ItoI up to a given order in the lattice spacing, a: 

Csym = C + aC^^^ + a'C^^^ + . . . , (2) 



® In fact, hairpin propagators will give larger exponential dependence than standard propagators because 
they are more chirally sensitive. 



For a related discussion, see Ref. 
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where contains higher dimensional operators of dimension 4 + n. Next one uses the 

method of spurion analysis to map the Symanzik action onto a chiral Lagrangian, in terms 
of pseudo-Goldstone mesons, which now incorporates the lattice spacing effects. This has 
been done in detail for a mixed GW-staggered theory in Ref. js^; here we only describe the 
results. 

The leading quark level Lagrangian is given by 

'iNf+2Ny 

C= Yl [^JP - ^Q\a Qb. (3) 

a,b=l 

where the quark fields are collected in the vectors 

Q^/=2 = ( ^ Jl,j2,j3,j4,^l,^2,^3,^ , (4) 

valence sea ghost 

QNf=2+l ^ ^^j^^ -^^ -^^ i^j^^ ^^^ ^ d, S)'^ (5) 

valence sea ghost 

for the two theories. There are 4 tastes for each flavor of sea quark, j, I, r.^° We work in the 
isospin limit in both the valence and sea sectors so the quark mass matrix in the 2+1 sea 
flavor theory is given by 

rriQ = diag(m„, rriu, m^, mj,mj, rrij, mj, rrij, mj,mj, rrij, rrir, rrir, rrir, rrir, mu, rriu, rris). (6) 

valence sea ghost 

The quark mass matrix in the two flavor theory is analogous but without strange valence, 
sea and ghost quark masses. The leading order mixed action Lagrangian, Eq. (jS)), has an 
approximate graded chiral symmetry, SU{ANf + N^\N^)l ® SU{ANf + N^\N^)ji, which is 
exact in the massless limit. In analogy to QCD, we assume that the vacuum spontaneously 
breaks this symmetry down to its vector subgroup, SU {AN f + N^\Ny)y , giving rise to {ANf + 
2Nyy — 1 pseudo-Goldstone mesons. These mesons are contained in the field 




^° Note that we use different labels for the valence and sea quarks than Ref. ■ Instead we use the "nuclear 
physics" labeling convention, which is consistent with Ref. [s^ . 

This is a "fake" symmetry of PQQCD. However, it gives the correct Ward identities and thus can be used 
to understand the symmetries and symmetry breaking of PQQCD [sH . 
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The matrices M and M contain bosonic mesons while x ci^nd contain fermionic mesons. 
Specifically, 





^ Vu 


7r+ .. 


■ 4>uj 


4>ui ■ 






n~ 


Vd ■■ 


■ <Pdj 


(pdl ■ 




M = 


4>ju 


^jd ■■ 




(pjl ■ 






4>iu 


4>id '■ 




Vi ■ 






\ '■ 

(<!>- 

H^uu 


4>ud ■ 






■J 


X = 




(Pdd ■ 


■ "^dj 


^dl ■ 





M 



^ V u 7r+ 

V • 
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(8) 



y : : : : : -.y 

In Eq. (jH} we only explicitly show the mesons needed in the two fiavor theory. The ellipses 
indicate mesons containing strange quarks in the 2+1 theory. The upper A^^^ x block 
of M contains the usual mesons composed of a valence quark and anti-quark. The fields 
composed of one valence quark and one sea anti-quark, such as (puj, are 1x4 matrices of 
fields where we have suppressed the taste index on the sea quarks. Likewise, the sea-sea 
mesons such as (pji are 4x4 matrix-fields. Under chiral transformations, S transforms as 



E — ^LER^ , L,Re SU{4Nf + N^\Ny)L,R. 



(9) 



In order to construct the chiral Lagrangian it is useful to first define a power-counting 
scheme. Continuum xPT is an expansion in powers of the pseudo-Goldstone meson momen- 



tum and mass squared 



41 



421: 



~ pI/ Al ml/Al 



(10) 



where oc rriQ and is the cutoff of xPT. In a mixed theory (or any theory which 
incorporates lattice spacing artifacts) one must also include the lattice spacing in the power 
counting. Both the chiral symmetry of the Ginsparg- Wilson valence quarks and the remnant 
U{1)a symmetry of the staggered sea quarks forbid operators of dimension five; therefore 
the leading lattice spacing correction for this mixed action theory arises at (9(a^). Moreover, 
current staggered lattice simulations indicate that taste-breaking effects (which are of (9(a^)) 
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are numerically of the same size as the lightest staggered meson mass Q • We therefore adopt 
the following power-counting scheme: 

~ Vll^l ~ "^q/Aqcd ~ a^A^cD • (11) 



The leading order (LO), 0{e^)^ Lagrangian is then given in Minkowski space by 38 1 

C = l^sii + ^str (^Sm^ + rngS^) - a" {Us + U's+ Uy) , (12) 

where we use the normalization / ~ 132 MeV and have already integrated out the taste 
singlet $0 field, which is proportional to str(<l>) lAs and lA'g are the well-known taste 
breaking potential arising from the staggered sea quarks 0, O]. The staggered potential 
only enters into our calculation through an additive shift to the sea-sea meson masses; we 
therefore do not write out its explicit form. The enhanced chiral properties of the mixed 
action theory are illustrated by the fact that only one new potential term arises at this order: 

Uv = -CMix str (rsSrsS"^) , (13) 

where 

r3=Vs-Vv = diag(- Jy, It ® Is, -Iv)- (14) 

The projectors, Vs and "Py, project onto the sea and valence-ghost sectors of the theory, /y 
and Is are the valence and sea fiavor identities, and It is the taste identity matrix. From 
this Lagrangian, one can compute the LO masses of the various pseudo-Goldstone mesons 
in Eq. (jH)). For mesons composed of only valence (ghost) quarks of fiavors a and 6, 

mlt^ = B{ma + m^). (15) 

This is identical to the continuum LO meson mass because the chiral properties of Ginsparg- 
Wilson quarks protect mesons composed of only valence (ghost) quarks from receiving mass 
corrections proportional to the lattice spacing. However, mesons composed of only sea 
quarks of fiavors Si and S2 and taste t, or mixed mesons with one valence (f ) and one sea 
(s) quark both receive lattice spacing mass shifts. Their LO masses are given by 

^ls,,t = B{m,, + + a2A(6), (16) 
ml^ = B{my + rris) + a^AMix- (17) 
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From now on we use tildes to indicate masses that include lattice spacing shifts. The 
only sea-sea mesons that enter tttt scattering to the order at which we are working are 
the taste-singlet mesons (this is because the valence-valence pions that are being scattered 
are tasteless), which are the heaviest; we therefore drop the taste label, t. The splittings 
between meson masses of different tastes have been determined numerically on the MILC 
configurations [^, so A(^/) should be considered an input rather than a fit parameter. The 
mixed mesons all receive the same shift given by 

Aiviix = — -p — , (18) 

which has yet to be determined numerically. 

After integrating out the $0 field, the two point correlation functions for the flavor-neutral 
states deviate from the simple single pole form. The momentum space propagator between 
two flavor neutral states is found to be at leading order 51 1 

- ^1 + (p2 - ml + ze) (p2 - ml + ze) UkQ\p^ " + ' 

where 

{+1 for a = valence or sea quarks 
(20) 
— 1 for a = ghost quarks . 

In Eq. (|19p . k runs over the flavor neutral states (0jj, 0^/, 0rr) and k' runs over the mass eigen- 

states of the sea sector. For tttt scattering, it will be useful to work with linear combinations 

of these rja flelds. In particular we form the linear combinations 

= ^{Vu - Vd) , ^ = + "^d) , (21) 

for which the propagators are 

QAp') = , (22) 



— + le 

i 



Speciflcally, 



i 2i{p'^ -rh^Mp^ -ml) 



UV") = -2 2 - '7^2 1^' for = 2, (24) 



p2 _ ^2 3 ^p2 _ fjiiy [p^ — ml) 
where = \{m^jj + 2ml). 



for Nf = 2 + 1, (25) 



11 



IV. CALCULATION OF THE 1 = 2 PION SCATTERING AMPLITUDE 



Our goal in this work is to calculate the / = 2 tttt scattering length in chiral perturbation 
theory for a partially quenched, mixed action theory with GW valence quarks and staggered 
sea quarks, in order to allow correct continuum and chiral extrapolation of mixed action 
lattice data. We begin, however, by reviewing the pion scattering amplitude in continuum 
SU (2) chiral perturbation theory. We next calculate the scattering amplitude in Nf = 2 
PQxPT and MA^PT, and finally in Nf = 2 + 1 PQxPT and MA^PT. When renormal- 
izing divergent 1-loop integrals, we use dimensional regularization and a modified minimal 
subtraction scheme where we consistently subtract all terms proportional to |4l| : 

2 

--7i5 + log47r + l, 

4 — ct 

where d is the number of space-time dimensions. The scattering amplitude can be related 
to the scattering length and other scattering parameters, as we discuss in Section El 



A. Continuum SU{2) 

The tree-level 1 = 2 pion scattering amplitude at threshold is well known to be 

. , Aiml 
lA = 







J n 



(26) 



It is corrected at 0(e^) by loop diagrams and also by tree level terms from the NLO (or 
Gasser-Leutwyler) chiral Lagrangian ^.^^ The diagrams that contribute at one loop or- 
der are shown in Figure P they lead to the following NLO expression for the scattering 
amplitude: 



iA 



Aim, 



■Pi=0 



P 



1 + 



81n' """" 



(27) 



where niuu is the tree-level expression given in Eq. (fTH|l and / is the LO pion decay constant 
which appears in Eq. 0. The coefficient l'^^ is a linear combination of low energy constants 
appearing in the Gasser-Leutwyler Lagrangian whose scale dependence exactly cancels the 
scale dependence of the logarithmic term. One can re-express the amplitude, however, in 
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The continuum tttt scattering amplitude is known to two-loops [43, |4J, |4 
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Pi 




(a) 



Pi 



P2 





(d) 



(b) 





(e) 



(c) 



FIG. 1: One- loop diagrams contributing to the vrvr scattering amplitude. Diagrams (a)-(c) are the 
S-, t-, and ti-channel diagrams, respectively, while diagram (e) represents wavefunction renormal- 
ization. 

terms of the physical pion mass and decay constant using the NLO formulae for and 
to find: 



J TT 



1 + 



mt, 



31n( ) + 



(28) 



where /,r7r is a different linear combination of low energy constants. The expression for /tttt can 
be found in Ref. [4^. We do not, however, include it here because we do not envision either 
using the known values of the Gasser-Leutwyler parameters in the the fit of the scattering 
length or using the fit to determine them. The simple expression 1)2811 has already been used 
in extrapolation of lattice data from mixed action simulations 29|, but it neglects lattice 
spacing effects from the staggered sea quarks which are known from other simulations to be 
of the same order as the leading order terms in the chiral expansion of some observables P| . 
We therefore proceed to calculate the scattering amplitude in a partially quenched, mixed 
action theory relevant to simulations. 



B. Mixed GW-Staggered Theory with two Sea Quarks 

The scattering amplitude in the partially quenched theory differs from the unquenched 
theory in three important respects. First, more mesons propagate in the loop diagrams. 
Second, some of the mesons have more complicated propagators due to hairpin diagrams at 

13 




FIG. 2: Example quark flow for a one- loop t-channel graph. This diagram illustrates the presence 
of meson loops composed of purely valence- valence mesons which are not canceled by valence-ghost 
loops. Different colors (shades of grey) represent different quark flavors. 



the quark level 



. Thi rd, there are additional terms in the NLO Lagrangian which 



IG 



38|. 



arise from partial quenching |52], and lattice spacing effects 

At the level of quark flow, there are diagrams such as Figure which route the valence 
quarks through the diagram in a way which has no ghostly counterpart. Consequently, 
the ghosts do not exactly cancel the valence quarks in loops. Of course, this is simply a 
reflection of the fact that the initial and final states — valence pions — are themselves 
not symmetric under the interchange of ghost and valence quarks, and therefore the graded 
symmetry between the valence and ghost pions has already been violated. This is well known 



in quenched and partially quenched heavy baryon xPT 







72 



73 . This fact also partly 

explains the success of quenched vrvr scattering in the 1 = 2 channel 2^ |^ ; quenching does 
not eliminate all loop graphs like it does in many other processes, and in particular, the 
s-channel diagram is not modified by (partial) quenching effects. As a consequence, it is 
necessary to compute all the graphs contributing to this process in order to determine the 
scattering amplitude. 

Quark level disconnected (hairpin) diagrams lead to higher ord er p oles in the propagator 
of any particle which has the quantum numbers of the vacuum In the isospin limit of 

the Nf = 2 partially quenched theory, conservation of isospin prevents the 7r° from suffering 
any hairpin effects. Hence only the f] acquires a disconnected propagator. Moreover, in the 
rriQ oo limit, the f] propagator (given for a general PQ theory in Eq. pHj) ) is given by the 
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FIG. 3: Example hairpin diagrams contributing to pion scattering. The propagator with a cross 
through it indicates the quark-disconnected piece of the fj propagator, Eq. H29|) . 



simple expression 



2 ~ 2 

— • 

— i- 



(29) 



(^p2 _ ^ 



2^2 ' 



where the parameter 



^lQ = rhl-ml (30) 

quantifies the partial quenching. (Recall that rhjj is the physical mass of a taste singlet 
sea-sea meson.) Notice that when Apg ^ the propagator also goes to zero; this is 
what we expect since, in the SU{2) theory, the only neutral propagating state is the 7r°. 
The propagator in Eq. can appear in loops, thereby producing new diagrams such as 
those in Fig. El^'^ After adding all such hairpin diagrams, one finds that the contribution of 
the fj to the amplitude is 

^A- - -J^^ (31) 

In addition to 1-loop contributions, the NLO scattering amplitude receives tree-level 
analytic contributions from operators of 0{e^) in the chiral Lagrangian. At this order, the 
mixed action Lagrangian contains the same 0(p^), 0{p'^mq), and 0{mg) operators as in 



We note that there are also similar contributions to the four particle vertex with a loop and to the mass 
correction. We do not show them, however, because they cancel against one another in the amplitude 
expressed in lattice-physical parameters, which we will show in the following pages. 

We note that this contribution does not vanish in the limit that with m^,- 7^ 0. Similar effects 

have been observed in quenched computations of pion scattering amplitudes |59l |60 ] . This non- vanishing 
contribution is the 1 = 2 remnant of the divergences that are known to occur in the / = amplitude at 
threshold. These divergences give rise to enhanced volume corrections to the 1 = amplitude with resp ect 
to the one-loop 1 = 2 amplitude and prevent the use of Liischer's formula. Moreover, it is known |48l l5Cl| 
that PQxPT is singular in the limit m„ with nonzero sea quark masses, so the behavior of the 
amplitude in this limit is meaningless. 
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the continuum partially quenched chiral Lagrangian, plus additional 0(a^), O{o?mq), and 
O{o?p'^) operators arising from discretization effects. We can now enumerate the generic 
forms of analytic contributions from these NLO operators. Because of the chiral symmetry 
of the GW valence sector, all tree-level contributions to the scattering length must vanish in 
the limit of vanishing valence quark mass.^^ Thus there are only three possible forms, each 
of which must be multiplied by an undetermined coefficient: m^^, m^^m^^-, and m^„a^. It 
may, at first, seem surprising that operators of 0{a^mq), which come from taste-symmetry 
breaking and contain projectors onto the sea sector, can contribute at tree-level to a purely 
valence quantity. Nevertheless, this turns out to be the case. These 0{a^mq) mixed action 
operators can be determined by first starting with the NLO staggered chiral Lagrangian 1^ , 
and then inserting a sea projector, Vs-, next to every taste matrix. One example of such an 



operator is str ^Sm^ j str (V s^^'^i^'^'^) + P-c 



, where, ^5 is the 75 matrix acting in taste- 



space and p.c. indicates parity-conjugate. This double-trace operator will contribute to the 
lattice pion mass, decay constant, and 4-point function at tree-level because one can place 
all of the valence pions inside the first supertrace, and the second supertrace containing the 
projector Vs will just reduce to the identity. 

Putting everything together, the total mixed action scattering amplitude to NLO is 



Aiml 



1 + 



m: 



TTl 

41n' 



+ 4 



%in 



(47r/)^ 



^PQ ^PQ 



6ml 



m: 



In 



+ 



(477/) 



(32) 



The first line of Eq. (j32p contains those terms which remain in the continuum and full QCD 
limit, Eq. (jTTj) . while the second line accounts for the effects of partial quenching and of 
the nonzero lattice spacing. Note that, for consistency with the 1-loop terms, we chose 
to re-express the analytic contribution proportional to the sea quark mass as m^^Apg. In 
Eq. (|32j) we have multiplied every contribution from diagrams which contain a sea quark 
loop by 1/4, thus making our expression applicable to lattice simulations in which the fourth 
root of the staggered sea quark determinant is taken. 

It is useful, however, to re-express the scattering amplitude in terms of the quantities 



As we discussed in the previous footnote, this condition need not hold for loop contributions to the 
scattering amplitude. 
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that one measures in a lattice simulation: and /,r- Throughout this paper, we will 
refer to these renormalized measured quantities as the lattice-physical pion mass and decay 
constant. Because we are working consistently to second order in chiral perturbation 
theory, we can equate the lattice-physical pion mass to the 1-loop chiral perturbation theory 
expression for the pion mass, and likewise for the lattice-physical decay constant. Thus, in 
terms of lattice-physical parameters, the mixed action J = 2 vrvr scattering amplitude is 



Aimi _ m. 



2 



31n(^^ ) -l + lUii) 



(47r/02g^4 



(33) 



where the first few terms are identical in form to the full QCD amplitude, Eq. (j28|) . This 
expression for the scattering amplitude is vastly simpler than the one in terms of the bare 
parameters. First, the hairpin contributions from all diagrams except those in Fig. |21 have 
exactly cancelled, removing the enhanced chiral logs and leaving the last term in Eq. ()33p as 
the only explicit modification arising from the partial quenching and discretization effects. 
Second, all contributions from mixed valence-sea mesons in loops have cancelled, thereby 
removing the new mixed action parameter, Cmix, completely. Third, all tree- level contri- 
butions proportional to the sea quark mass have also cancelled from this expression. And 
finally, most striking is the fact that an explicit computation of the 0{a'^mq) contributions 
to the amplitude arising from the NLO mixed action Lagrangian show that these effects 
exactly cancel when the amplitude is expressed in lattice-physical parameters. This result 
will be discussed in detail in Ref. 7^. Thus to reiterate, the only partial quenching and lat- 
tice spacing dependence in the amplitude comes from the hairpin diagrams of Fig. El which 
produce contributions proportional to Apg = (m^^- + a^A(^/) — m^)^, where m'jj + a^A(^i) 
is the mass-squared of the taste-singlet sea-sea meson. Moreover, we presume that anyone 
performing a mixed action lattice simulation will separately measure the taste-singlet sea- 
sea meson mass and use it as an input to fits of other quantities such as the vrvr scattering 
length. Thus we do not consider it to be an undetermined parameter. 



Notice that once the lattice spacing a has been determined, the lattice-physical pion mass can be unam- 
biguously determined by measuring the exponential decay of a pion-pion correlator. We assume that the 
lattice spacing a has been determined, for example, by studying the heavy quark potential or quarkonium 
spectrum. 

Another consequence of the exact cancelation of the loops with mixed valence-sea quarks is that one does 
not have to implement the "fourth-root trick" through this order. 
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It is now clear that one should fit vrvr scattering lattice data in terms of the lattice- 
physical pion mass and decay constant rather than in terms of the LO pion mass and LO 
decay constant. By doing this, one eliminates three undetermined fit parameters: Cmix, I'pq, 
and /' 2 , as well as the enhanced chiral logs. 



C. Mixed GW-Staggered Theory with 2 + 1 Sea Quarks 

The 2 + 1 fiavor theory has three additional quarks - the strange valence and ghost 
and strange sea quarks - which can lead to new contributions to the scattering amplitude. 
Because we only consider the scattering of valence pions, however, strange valence quarks 
cannot appear in this process. Thus all new contributions to the scattering amplitude 
necessarily come only from the sea strange quark, r. Because the r quark is heavier than 
the other sea quarks there is SU (3) symmetry breaking in the sea. This symmetry breaking 
only affects the pion scattering amplitude, expressed in lattice-physical quantities, through 
the graphs with internal f] propagators because the masses of the mixed valence-sea mesons 
cancel in the final amplitude as they did in the earlier two fiavor case. In addition, the 
only signature of partial quenching in the amplitude comes from these same diagrams. It is 
therefore worthwhile to investigate the physics of the neutral meson propagators further. 

There are more hairpin graphs in the 2 + 1 fiavor theory since the rjg may propagate as 
well as the rju and the rjd. Because these mesons mix with one another, the fiavor basis is 
not the most convenient basis for the computation. Rather, a useful basis of states is vr", 
V = iVu + Vd) I and rig. Since we work in the isospin limit, the 7r° cannot mix with or //<,; 
in addition, there is no vertex between the ri^ and vr+vr^ at this order, so we never encounter 
a propagating r^^. Thus all the PQ effects are absorbed into the f] propagator, which is given 
by 

^ — 3 {p^ — m'^y [p^ — m^) 

In SU{'i) chiral perturbation theory, the neutral mesons are the vr*^ and the 773. Therefore, 
in the PQ theory, we know that there will be a contribution from the f] graphs that does not 
result from partial quenching or SU (3) symmetry breaking. Therefore the extra PQ graphs 
arising from the internal fj fields must not vanish in the Apg ^ limit, in contrast to the 
two fiavor case of Eq. (|^ . 
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To make this clear, we can re-express the propagator of Eq. ()34|) in terms of Apg as 



m 



2\2 



1 

+ 



A2 



3j92 



rat 



(35) 



This propagator has a single pole which is independent of Apg, as well as higher order poles 
that are at least quadratic in Apg. It is interesting to consider the large limit of this 
propagator. In this limit, ^ ^Brrir is also large. For momenta that are small compared 
to m^, the second term of this equation goes to zero in the large limit, and the 2 + 1 
flavor propagator reduces to the 2 flavor propagator, Eq. (j^ . as expected. 

While the above expression clarifles the Apg dependence of the propagator and the large 
iTLr limit, it obscures the SU{'i)sea limit. An equivalent form of the propagator is 



^PQ 
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ni: 
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^PQ 
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(36) 



where the quantity Aj 



rrij^ parametrizes the SU{'i)sea breaking. When A3 = 



this propagator is similar in form to the corresponding 2 flavor propagator, Eq. (j^^ . but it 
has an additional single pole due to the extra neutral meson in the SU (3) theory. 

Having considered the new physics of the hairpin propagator, we can now calculate the 
scattering amplitude. For our purposes here, it is most convenient to express the total 
J = 2 TTvr scattering amplitude in terms of Apg. Just as in the 2-flavor computation, 
the NLO analytic contributions due to partial quenching and flnite lattice spacing effects 
exactly cancel when the amplitude is expressed in lattice-physical parameters. All sea quark 
mass and lattice spacing dependence comes from the hairpin diagrams, which produce terms 
proportional to powers of Apg with known coefficients. The amplitude is 
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19 



where Apg = rrijj + a^A(^/) — ml and 

•^1^^) = - 9(3.^1)2 - 1) - (3a^ + 2) ln(x)] , (38a) 

^2(x) = g^^^^^3 [(x - l)(x + 3) - (3x + 1) ln(x)] , (38b) 

= 9(^i^)4 [(^ - - 7a; - 12) + 2{7x + 2) ln(x)] , (38c) 

J^i{x) = - ^ [(x - l)(x2 - 8x - 17) + 6(3x + 1) ln(x)] . (38d) 

The functions jFj have the property that J-'i{x) — ^> in the hmit that x —>■ 00. Therefore, when 
the strange sea quark mass is very large, i.e. m^/ml ^ 1, the 2 + 1 flavor amplitude reduces 
to the 2 flavor amplitude, Eq. (jHHj) . with the exception of terms that can be absorbed into 
the analytic terms. The low energy constants have a scale dependence which exactly cancels 
the scale dependence in the logarithms. The coefficient /tttt is the same linear combination 
of Gasser-Leutwyler coefficients that appear in the SUCi) scattering amplitude expressed in 
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terms of the physical pion mass and decay constant 

Because the functions jFj depend logarithmically on x, the 2+1 flavor scattering amplitude 
features enhanced chiral logarithms |4S| that are absent from the 2 flavor amplitude. This is 
a useful observation, as we will now explain. Because there is a strange quark in nature and 
its mass is less than the QCD scale, Aqc/), lattice simulations must use 2 + 1 quark flavors. 
It is often practical to flx the strange quark mass at a constant value near its physical value 
in these simulations. This circumstance is helpful because, just as SU{2) chiral perturbation 
theory is useful to describe nature at scales smaller than the strange quark mass, the 2 flavor 
amplitude given in Eq. ()33|1 can be used to extrapolate 2 + 1 flavor lattice data at energy 
scales smaller than the strange sea quark mass used in the simulation (provided, of course, 
there are no strange valence quarks) [2^. This is valid because, at energy scales smaller than 
the strange quark mass (or actually twice the strange quark mass, since the purely pionic 
systems have no valence strange quarks), one can integrate out the strange quark. This is 
not an approximation, because all of the effects of the strange quark are absorbed into a 
renormalization of the parameters of the chiral Lagrangian. Moreover, since the 2 flavor 
amplitude does not exhibit enhanced chiral logarithms, signatures of partial quenching can 
be reduced by extrapolating lattice data with the 2 flavor, rather than the 2 + 1 flavor, 
expression. We note that in this case the effects of the strange quark are absorbed in the 
coefficients of the analytic terms appearing in Eq. and thus they are not constant, but 
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rather depend logarithmically upon the strange sea quark mass. 



V. 1 = 2 PION SCATTERING LENGTH RESULTS 



In this section we present our results for the s-wave / = 2 tttt scattering length in the two 
theories most relevant to current mixed action lattice simulations: those with GW valence 
quarks and either Nf = 2 or Nf = 2 + 1 staggered sea quarks. We only present results for 
the scattering length expressed in lattice-physical parameters. The s-wave scattering length 
is trivially related to the full scattering amplitude at threshold by an overall prefactor: 

327rm^ 

A. Scattering Length with 2 Sea Quarks 

The / = 2 TTTT s-wave scattering length in a MA^PT theory with 2 sea quarks is given by 



(39) 



87r/2 {ATiUy 



31n( ) -l + /^^(^) 



ml At 



PQ 



(40) 



where Apg = m'jj + a'^A{^j) — m^. The first two terms are the result one obtains in SU{2) 
xPT and the last term is the only new effect arising from the partial quenching and mixed 
action. All other possible partial quenching terms, enhanced chiral logs and additional linear 
combinations of the 0{p^) Gasser-Leutwyler coefficients, exactly cancel when the scattering 
length is expressed in terms of lattice-physical parameters. And, most strikingly, the pion 
mass, decay constant and the 4-point function all receive 0{a^mq) corrections from the 
lattice, but they exactly cancel in the scattering length expressed in terms of the lattice 
-physical parameters 2^. It is remarkable that the only artifact of the nonzero lattice 



spacing, m^-- + a'^Aj, can be separately determined simply by measuring the exponential 
fall-off of the taste-singlet sea-sea meson 2-point function. Thus there are no undetermined 
fit parameters in the mixed action scattering length expression from either partial quenching 
or lattice discretization effects; there is only the unknown continuum coefficient, /tttt- 

One can trivially deduce the continuum PQ scattering length from Eq. ()40|1 : simply let 
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0, reducing m 



m 



2Bmj in Apg, resulting in 
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(41) 



B. Scattering Length with 2+1 Sea Quarks 



The I = 2 TTTT s-wave scattering length in a MA^PT theory with 2+1 sea quarks is given 



by 
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(42) 



where the functions JFj are defined in Eq. (j38|) . As in the 2-flavor MA^PT expression, 
Eq. (j4Up . the only undetermined parameter is the linear combination of Gasser-Leutwyler 
coefficients, /tttt, which also appears in the continuum xPT expression. 

We note as an aside that this suppression of lattice spacing counterterms is in contrast 
to the larger number of terms that one would need in order to correctly fit data from 
simulations with Wilson valence quarks on Wilson sea quarks. Because the Wilson action 
breaks chiral symmetry at 0{a), even for massless quarks, there will be terms proportional 
to all powers of the lattice spacing in the expression for the scattering length in Wilson 
XPT ll^ 



761, 



77|. Moreover, such lattice spacing corrections begin at 0{a), rather than 



0(0^). If one uses 0{a) improved Wilson quarks, then the leading discretization effects 
are of O(a^), as for staggered quarks; however, this does not remove the additional chiral 
symmetry-breaking operators. Another practical issue is whether or not one can perform 
simulations with Wilson sea quarks that are light enough to be in the chiral regime. 



VI. DISCUSSION 



Considerable progress has recent 



scattering in the 1 = 2 channel 
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y been made in fully dynamical simulations of pion 
" \. We have considered 1 = 2 scattering of pions 
composed of Ginsparg- Wilson quarks on a staggered sea. We have calculated the scattering 
length in both this mixed action theory and in continuum PQ^PT for theories with either 



22 



2 or 2 + 1 dynamical quarks. These expressions are necessary for the correct continuum and 
chiral extrapolation of PQ and mixed action lattice data to the physical pion mass. 

Our formulae, Eqs. (|4Up . (j42j) . not only provide the form for the mixed action scattering 
length, but also contain two predictions relevant to the recent work of Ref. ^qI. Beane 
et. al. calculated the 1 = 2 s— wave tttt scattering length using domain wall valence quarks 
and staggered sea quarks, but used the continuum xPT expression to extrapolate to the 
physical quark masses. In Figure 2 of Ref. [2^, which plots m-^af^ versus rriT,/ f^,, the fit of 
the xPT expression to the lattice data overshoots the lightest pion mass point but fits the 
heavier two points quite well. This is interesting because Eq. pn|l predicts a known, positive 
shift to m^a^^ of size ApQ/(768/^7r^). Accounting for this positive shift is equivalent to 
lowering the entire curve, and could therefore move the fit such that it goes between the 



data points. This turns out, however, not to be the case. In Ref. j29L the valence and sea 
quark masses are tuned to be equal, so Apg = a^A/ ~ (446 MeV)^ 5|. Despite the large 
value of ApQ, the predicted shift is insignificant, being an order of magnitude less than the 
statistical error. In Table HI we collect the predicted shifts to mT^af^ at the three pion masses 
used in Ref. We also list the magnitude of the ratio of these predicted shifts to the 

leading contribution to the scattering length, which turn out to be small, lending confidence 
to the power counting we have used, Eq. (fTT|) . The other more important prediction is 
that there are no unknown corrections to the xPT formula for the scattering length arising 
from lattice spacing corrections or partial quenching through the order 0{m'^), 0{a^mq) and 
0{a^). Therefore, to within statistical and systematic errors, the continuum xPT expression 
used by Beane et. al. to fit their numerical tttt scattering data [3| receives no corrections 
through the 1-loop level. 

The central result of this paper is that the appropriate way to extrapolate lattice tttt 
scattering data is in terms of the lattice-physical pion mass and decay constant rather than in 
terms of the LO parameters which appear in the chiral Lagrangian. When expressed in terms 
of the LO parameters, the scattering length depends upon 4 undetermined parameters, Z^^, 
I'pQ, l'^2 and CMix- In contrast, the scattering length expressed in terms of the lattice-physical 
parameters depends upon only one unknown parameter, /tttt; the same linear combination of 
Gasser-Leutwyler coefficients which contributes to the scattering length in continuum xPT. 
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TABLE I: Predicted shifts to the scattering length computed in Ref. 29| arising from finite lattice 
spacing effects in the mixed action theory. The first two rows show the approximate values of m,r 



and while the third shows mT^a"2^ plus the statistical error calculated in 
we give the predicted shifts in the scattering length (times ttt-tt) and, in the fifth row, we give the 
ratio of the predicted shift to the leading order contribution to the scattering length. 



29|. In the fourth row. 



(MeV) 
U (MeV) 

(0) 


294 348 484 
145 149 158 
-0.212 ± 0.024 -0.222 ± 0.014 -0.38 ± 0.03 


A* 

76877^/^ 


0.00374 0.00336 0.00266 
0.0229 0.0155 0.00711 
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